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^ Abstract 

I discuss radiative corrections due to the emission of photons from quarks 
which contribute to deep inelastic lepton-nucleon scattering as well as to 
pp collisions at high energies. These corrections are dominated by quark- 
mass singularities which have to be absorbed into the parton distribution 
functions. Observable effects appear as a modification of the Q 2 dependence 
of the distribution functions. Numerical results indicate, however, that these 
QED corrections are negligible except at extremely large Q 2 and large x. 
Therefore it is safe to neglect the single and multiple photon effects in pp 
scattering at LHC energies. 
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1 Introduction 



High-energy scattering experiments with charged particles require the inclusion of 
electromagnetic radiative corrections due to the virtual and real emission of photons. 
These QED corrections are known to a high precision as far as they are related to ra- 
diation from leptons, as e.g. in the case of initial state radiation in e + e~ annihilation 
(see, e.g. [1]) or in the case of leptonic corrections in deep inelastic lepton-nucleon 
scattering [2]. In processes involving charged hadrons in the initial or final state 
there are QED corrections related to photon emission from hadrons or their con- 
stituents. Photons originating from hadron decays or emerging during hadronization 
are usually modeled in Monte Carlo programs simulating the hadronic final state 
and are not considered in this paper. In the following I want to discuss the effect 
of photon radiation from quarks entering or leaving the underlying hard scattering 
process. I will show that the corresponding corrections are negligible, even at very 
high energies. 

In e + e~ annihilation into hadrons, QED corrections of the type considered here 
appear as photon radiation from the final state quarks. Final state radiation does not 
pose a severe problem since as a consequence of the KLN theorem the corresponding 
corrections are small. The situation is different in deep inelastic scattering or in 
hadron collisions. In this case there are QED corrections related to the emission of 
photons from quarks in the initial state. These corrections contain mass singularities 
related to the initial state partons. Using dimensional regularization — as usual in 
QCD calculations — these mass singularities appear as poles in 1/(D — 4). If one 
regulates the mass singularities with the help of finite quark masses — as wide-spread 
in QED calculations — one finds terms proportional to lnm 2 . 

In an early calculation of radiative corrections to charged current neutrino scat- 
tering by Kiskis [3], the author was worried by this fact, since it makes a large 
difference whether one uses m q = 350 MeV (a typical constituent quark mass) lead- 
ing to In /i 2 /m 2 ~ 10.9 for /i 2 = M^, or m q = 5 MeV (a typical value for the 
current quark mass) leading to In M^/m 2 ~ 19.4. Later, De Puijula, Petronzio and 
Savoy-Navarro [4] have argued that the unphysical dependence on quark masses can 
be absorbed by a redefinition of parton distribution functions. This redefinition, 
or renormalization of parton distribution functions, is well-known in the calculation 
of QCD radiative corrections where in complete analogy to photon radiation the 
emission of gluons leads to mass singularities as well. 

The essential property of mass singularities is that they factorize [5], i.e. they 
can be written as a convolution of the parton-level Born cross-section o v ^ r or ° n n with a 
kernel A"(^,/i 2 ) describing the effect of radiation and containing the mass-singular 
terms. Explicitly, for first-order corrections and omitting summation over parton 
types, the mass-singular part has the form 

f fy/fO^^VC^P). (1) 

Z7T Jx Q 

Here /(£) is a parton distribution function, p the momentum of the incoming hadron 
and £ a dimensionless variable characterizing the amount of energy available for the 
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hard scattering process after radiation. The minimum of £, ^ m in = x, is determined 
by the kinematics of the process. The energy scale ji (factorization scale) remains 
arbitrary unless the non-singular contributions are specified. Choosing ji to be a 
typical mass scale of the process, it is usually possible to avoid the appearance of 
large non-singular corrections. The lowest-order cross section is itself a convolution 
of parton cross sections with parton distribution functions, 



Adding these contributions, one can absorb the mass-singular terms by a redefinition 
of the parton distribution functions: 



—5({ - x ) — 



a 



8(l-u) + -^-K{tlu,n 2 ) 
Ztt 



(3) 



s(i-x)r-(i^)az:° n n (iP) 



with renormalized distribution functions 

f 1 du 

u 



rren{£ 2 



a 



8(l-u) + —K{tlu,n 2 ) 
Ztt 
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(4) 



The bare parton distribution functions /(£) are, in fact, not measurable. It is 
rather the renormalized distributions f ren that have to be identified with measured 
(i.e. finite) quantities. The appearance of mass singularities in unphysical quanti- 
ties is an artifact of the perturbative treatment. By the redefinition Eq. (4), the 
mass singularities disappear from the observable cross section and the renormalized 
distribution functions become dependent on the factorization scale /i 2 . The /i 2 de- 
pendence is controlled by the well-known Gribov-Lipatov-Altarelli-Parisi (GLAP) 
equations [6]. In deep inelastic scattering they are equivalent to the renormaliza- 
tion group equations for the Wilson coefficients and express the fact that observable 
effects are independent of the energy scale at which the distribution functions are 
renormalized. The solution of the GLAP equations corresponds to the resummation 
of all powers n of the leading logarithms a™(ln/i 2 ) n . 

Since mass singularities are universal, i.e. independent of the process under 
consideration, the definition of renormalized parton distributions is also universal. 
Therefore it is possible to discuss the bulk of radiative corrections in terms of par- 
ton distribution functions. This will be true if there is only one large scale in the 
process (e.g. in inclusive deep inelastic scattering ep — > eX with x not small). Then 
radiative corrections which are not mass-singular cannot contain large logarithms. 

The above prescription for the treatment of mass singularities applies to both 
QCD and QED corrections. Taking into account QED corrections, the renormal- 
ization of parton distributions has to include terms due to the emission of photons 
from quarks, in addition to those due to the emission of gluons. Then, also the 
GLAP equations are modified by an additional term of the order of the electro- 
magnetic fine-structure constant, a e . Apart from small non-singular contributions, 
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the resulting modified scale dependence of parton distribution functions is the only 
observable effect of QED corrections in high-energy scattering of hadrons. 

The modification of the GLAP evolution equations by QED corrections has also 
been discussed in [7]. In this work, I will present and discuss numerical results 
for their solutions including terms of order 0(a e (a s ln/i 2 ) n ) for arbitrary n relevant 
for present and future experiments like at HERA or at LHC. Moreover, I will give 
a simple prescription to approximately take into account the QED corrections of 
valence parton distributions. 

The prescription to be described below applies to completely inclusive measure- 
ments of any process hih 2 — >■ X } i.e. where emitted photons are not restricted by 
energy or angle cuts to a specific phase-space region. Experimentally this means 
that no attempt is made to observe emitted photons. It is a more complicated task 
to derive cross sections for the production of isolated hard photons (see for example 
[8] and references therein). In this case, mass-singular contributions can be avoided 
by imposing isolation cuts. For a measurement of direct photon production without 
isolation cuts one has to absorb part of the mass-singular contributions into pho- 
ton fragmentation functions describing the non-perturbative emission of collinear 
photons from partons. 

2 Formalism 

After having absorbed mass-singular terms into the parton distribution functions, 
qj(x } fi 2 ) for quarks with flavor / and G(x } fi 2 ) for gluons, the resulting dependence 
on the energy scale /i 2 at which the process probes the parton content of the hadrons 
is described by the GLAP evolution equations [6]. Using the scale variable 



t = ln/i 2 /A 2 , 



(5) 



they read 



d 
di 



q f (x,t) 




/•l dz r -i 

/ —\P q/q (z,t)q f (x/z,t) + P q/ cj(z,t)G(x/z,t)\, (6) 

Jx Z 



d 

di 



G(x,t) 




rl dz 

/ — Y. P G/ q (^thf(x/^t) + PG/G(z,t)G(x/z,t) . (7) 
Jx Z , 



To leading logarithmic accuracy (LLA), the splitting kernels P^j are scale inde- 
pendent and given by (see e.g. [10]) 
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P G/q (z) = C F l + (1 Z)2 ; 



P G/G (z) = 2N C 
1 



+ 



P q /rAz) 



z 2 + (l-z) 2 



In the following, I restrict myself to the LLA and omit the argument t in Pi/j. 
The running strong fine-structure constant is given by 



1 7 33 - 2N f 



(9) 



Note that the "+" distribution appearing in Eq. (8) is defined for the interval from 
to 1: 

f dzD + (z)f(z) = f dzD(z) (f(z) - f(l)) ; (10) 
JO Jo 

thus, if the integral is restricted to the range x < z < 1, an additional term 

/(l)ln(l — x) has to be taken into account. 

The inclusion of QED corrections modifies the evolution equation for the charged 

parton distributions by an additional term: 



d cts(t) Z" 1 dz r 

= I — \ P q/q( Z , t ) < lf( X / Z , t ) + Pq/G(z,t)G(x/z,t) 



dt 



2tt Jx z 

2lT Jx Z 



:iv, 



P q/q{ Z ^hf{ X / Z ^) 



where 



[l- z ) + 2 { ) 



q/q- 



'12) 



Here the running electromagnetic fine-structure constant appears, which is given by 

a(0) 



a e (t) 



1 



(13) 



e)(t-t mf )9(t-t mf ) 

with t m; = ln(mj/A 2 ) and a(0) = 1/137.036 .... ej are the fermion charges in units 
of the positron charge and the fermion-mass thresholds have been approximated by 
the step function 6{t — t m; ). 

It convenient to use the following combinations of parton distribution functions: 

^2 ( u ( x ) + u(x)) , 

gen 

(d(x) + d(i 

gen 

U(x) + D(x), 



U(x) 
D(x) 
E(x) 
A(x) 



(14) 



U(x) - D(x), 
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where J2 g en means summation over the generations. Then, denoting the convolution 
of splitting kernels with distribution functions by the symbol ® and the derivative 
with respect to the scale variable t by a dot, the evolution equations can be written 
as: 




In contrast to the pure QCD evolution, the non-singlet combination of parton distri- 
butions, A(x), does not decouple from the singlet sector due to the different charges 
of up- and down-type quarks. For our purpose it is convenient to separate the 
solution of the equations without the QED term, i.e. to write 

£ = E° + cr, A = A + 8, G = G° + g (16) 

where E°, A , and G° obey the evolution equations Eq. (15) with a e = 0. cr, 8, 
and g are corrections of relative order 0(a e ). Neglecting terms of order 0(a 2 e ), one 
obtains equations for the QED contributions to the parton distributions: 




Note that although the gluons do not couple directly to photons, their distribution 
function is modified through higher-order contributions induced by their coupling 
to the quarks. The correction g is formally of order 0(a e a s ). In Eqs. (17) the full 
a s evolution of the QED contributions is kept. 

3 Numerical Results 

The equations Eqs. (17) lend themselves directly to an iterative numerical solu- 
tion, given the QCD evolved solutions E°(x,t) and A°(x,t). The use of Eqs. (15) 
would require initial conditions for qj(x } t ) and G(x } t ) at some reference scale 
t = ln(Qo/A 2 ). High precision and a stable numerical algorithm is needed if the 
small corrections of 0(a e ) were to be determined directly from Eqs. (15) instead 
of Eqs. (17). Therefore it is preferable to use Eqs. (17). I checked my algorithm 
by comparing the solutions from these two equations. As an additional check, the 
algorithm was used to solve Eqs. (15) with a e = and initial conditions qj(x } t ) 
taken from one of the commonly used parametrizations of parton distribution func- 
tions, as provided for example by the program library PAKPDF [11]. The resulting 
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Q 2 dependence of the pure QCD evolution was then compared with the Q 2 depen- 
dence of the corresponding parametrization. In view of the fact that the available 
parametrizations are only approximate solutions of the GLAP equations having their 
own limited precision, the agreement between these different methods was remark- 
ably good. 

The solution of the differential equations Eqs. (15) or (17) requires initial condi- 
tions, i.e. the knowledge of parton distribution functions at a reference scale ji 2 = Qq. 
The input distribution functions qj(x,Q^), G(x } Qq) have to be taken from experi- 
mental data which (with only a few exceptions) have not been corrected for radiative 
corrections due to photon emission from quarks. Consequently, the QED corrections 
cr, 8 } and g in Eqs. (17) have to vanish at the reference scale. 

In Figs. 1, 2, and 3, I show results for the corrections to the distribution functions 
U(x } Q 2 ) } D(x } Q 2 ) } G(x } Q 2 ) and the structure function 



measured in deep inelastic electron-proton scattering at HERA. The figures show 
the QED corrections in per cent relative to the distribution functions obtained from 
the GLAP equations without the QED term. The input distributions were taken 
from [12] (set 1.1) since those parametrizations are fairly simple and the numerical 
solution of Eqs. (17) requires less computer time than for other parton distributions. 
The figures show small, negative corrections at the per-mille level for all values of 
x and Q 2 relevant in forthcoming experiments. Only at large x ~ 0.5 and large 
Q 2 ~ 10 3 GeV 2 the corrections reach the magnitude of one per cent. 

The increase of corrections for x — > 1 is due to the ln(l — x) terms appearing 
in the evaluation of the "+" distributions. Although large relative corrections are 
reached in this limit, they are of no practical relevance, since parton distributions 
and cross sections decrease as x — > 1 and to reach an experimental accuracy of a 
few per cent at large x will be very unlikely. 

The largest corrections are obtained for up-type quarks due to the larger charge 
factor 4/9 as compared to 1/9 for down-type quarks. Therefore the correction to 
F% is close to that for U(x). The gluon distribution, being of order 0(a e a s ) } is 
corrected by less than 0.1% up to values of x of about 0.2. 

Vice versa, the effect of g being non-zero on the quark distributions is of the 
same size. Fixing g at g(x) = would lead to different results, particularly at small 
x, though at the same level of a few per mille and with the same cross features. 
At the level of the observed QED corrections cr, 8, and g, their QCD evolution is 
non-negligible. 

The corrections vanish for Q 2 — > Qq since I assumed that the input distributions 
qj(x } t ) and G(x } t ) have been extracted from experiment at the reference scale Qq 
without subtracting quarkonic QED corrections. 

The asymptotic behavior for x — > can be checked analytically. It is well-known 
that a singular behavior of distributions oc x~ a for x — > remains unchanged by the 
GLAP equations if a > 1. Thus the 0(a e ) corrected distributions have the same 
power behavior as the uncorrected ones, the ratio consequently reaching a constant 
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value for x — > 0. The value of this constant depends on the power a and other details 
of the distribution functions. The valence parts of U(x) and D(x), however, which 
vanish at x = 0, receive positive corrections at small x, thus producing the well- 
known physical picture: radiation of gluons as well as of photons leads to a depletion 
at large x and an enhancement at small x, i.e. partons are shifted to smaller x. 

The corrections strongly depend on the input distribution functions; differences 
at the per-mille level are found. But still these differences are irrelevant when com- 
pared with the expected experimental precision of structure-function measurements. 

For valence distributions there is a simple prescription to include the effect of 
QED radiation. Using as variable the evolution length (in leading order) 

Ut)= / = —]*-, (19) 

Jto Zir ZTTOq to 

the GLAP equation for valence quark distributions has the simple form 



Then, the substitution 



d dz 

qval(x,is)= Pq/q(z,€s)qval(x/z,£ s ). (20) 



\ u ,(a s (f) t e) a e (t'Y 



Jtn \ ZlT 



t \ Z7T C'f ZlT J ^-p 



2C F J2 f e 2 f a e {t o y 

automatically takes account of QED radiative effects. This means, the evolution 
length is increased by QED radiation. Given the solution q va i(x,t) of Eq. (20), one 
obtains the solution of the GLAP equation including QED effects by 

q va i(x,t) q va i(x,i) (22) 

with 

t = t[ ^4 . 23 

This simple prescription is sufficient for including quarkonic QED corrections into 
the structure functions at large x where valence contributions dominate. 



4 Concluding Remarks 

I have shown that radiation of photons from quarks leads to negligible corrections 
of cross sections for high-energy scattering processes involving hadrons. The leading 
effect is a modification of the scale dependence of parton distribution functions 
described by a modification of the GLAP evolution equations. These equations 
include the effect of multiple photon emission as far as it is enhanced by a logarithm 
of the energy scale of the process. Non-leading corrections have not been calculated, 
but they are of the order of 0(a e /ir) ~ 0.2% and therefore also negligible. 
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A calculation of QED corrections using definite values for quark masses [13] 
and not factorizing the corresponding mass-singular contributions is misleading and 
usually leads to a gross overestimation of the effect. 

I did not discuss the question of how to model the production of soft or collinear 
photons in the final state of a hard scattering process. This should be viewed as a 
part of the non-perturbative hadronization process. Commonly used Monte Carlo 
programs apply a two-step procedure where in a first step a parton cascade is devel- 
oped. This cascade obeys the GLAP evolution and may allow for photon emission in 
addition to gluon radiation. The termination of the parton cascade requires a cut-off 
which could be implemented in a reasonable way by using effective quark masses. 
Usually those quark masses are identified with the constituent quark masses. It is 
important to understand, however, that this prescription does not correct the hard 
scattering cross section. The purpose of these Monte Carlo programs is to provide 
a model for the distribution of energy and momentum among final state hadrons. 
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Figure 1: Q 2 dependence of the QED corrections (in per cent, see text) to parton dis- 
tributions and the structure function F% for deep inelastic lepton-proton scattering 
at x = 0.1. Input parton distributions were taken from [12]. 
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Figure 2: x dependence of the QED corrections (in per cent, see text) to parton dis- 
tributions and the structure function F% for deep inelastic lepton-proton scattering 
at Q 2 = 6 X 10 3 GeV 2 . Input parton distributions were taken from [12]. 
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Figure 3: Q 2 dependence of the QED corrections (in per cent, see text) to the struc- 
ture function F% for deep inelastic lepton-proton scattering at x = 0.001, x = 0.1 
and x = 0.505. Input parton distributions were taken from [12]. 
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